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Large Multiplicity Gluon Production in QCD 



We compute the approximate cross section a n for producing re — 2 resolved 
gluons in a gluon-gluon collision, using the Parke- Taylor formula regularized in a 
Lorentz invariant manner. We find, in double leading logarithm approximation, 
that 



where ^/s cu t is the minimum invariant mass for a resolved gluon pair. There is no 
factor of l/(n — 2)! multiplying the expression. We present additional numerical 
results, and comment on their implications for perturbative calculations of re-jet 
cross sections at colliders. 
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Abstract 




Introduction 



There is little doubt that QCD, the SU(3) gauge theory describing the interactions 
of quarks and gluons, is the underlying theory of strong interactions. Various experi- 
mental results from both electron-positron machines |]J and hadron colliders M are in 
good agreement with QCD predictions. Use of sophisticated Monte Carlo generators 
based on coherent branching processes which can take into account both initial and 
final state radiation || are able to reproduce the CDF data Q up to six jets. 

Multiple QCD jets will, of course, constitute important background for new physics 
discovery at the future SSC/LHC colliders. Thus, there is practical importance in be- 
ing able to calculate differential cross sections for n-jet processes in QCD at collider 
energies, with cuts appropriate to experiments and detectors of interest. At a more 
theoretical level, recent work on the behavior of multiparticle amplitudes in field the- 
ories of massive bosons [0, ||, [7|, || || has revealed non-perturbative behavior for tree 
level cross sections near threshold. It is an interesting question whether the most es- 
tablished bosonic field theory (QCD) maintains a perturbative behavior at tree level 
for multiple jet production. 

In this Letter, we will focus on the behavior of the exclusive cross section for the pro- 
cess gg — > (n — 2) g All final state partons will be separated in a Lorentz-invariant 
way, by a given fixed minimum amount {Pi + Pj) 2 > s cut . We will dispense with regions 
of phase space where our cuts are not obeyed, so that we do not merge unresolved 
gluons into jets; it is in this sense that our result will be an exclusive cross section. We 
will concentrate on the limiting case of the production of a large number of final state 
gluons. Our results will depend on the ratio A = s cut /s, where s is the parton-parton 
center-of-mass (cm.) energy squared (usually denoted by s). Different values of A 
discriminate whether or not the produced jets constitute important backgrounds for 
new physics; it is also possible that multi-minijet jet cross sections (small A) may be 
the most important QCD contribution to the total hadronic cross section. A typical 
value of A for interesting jets at the SSC (pxi > 50 GeV, 9^ > 30°, t/s ~ 4 TeV) is 
A w lO -5 . 

We will also be concerned with the validity of perturbation theory in the approxi- 
mate tree-level calculation we will perform. Naively, each time there is an extra parton 
in the final state one should pay a price of a coupling constant. However, it is possible 
that in calculating certain quantities the expansion parameter becomes effectively the 
coupling constant multiplied by a number that may be large. In this case, care must be 
taken either by summing up these large contributions to all orders or, more modestly, 
performing perturbative calculations only when the effective expansion parameter is 
small. 

Exact multi-parton QCD amplitudes may be generated through use of the Berends- 
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Giele recursion relations [O]; amplitudes for n = 8 (six final state gluons) have been 



given explicitly ||12|| . However, the complexity of these amplitudes makes their useful- 
ness rather limited by the computer time required for their evaluation. Since it is the 
aim of this work to examine cross sections for large numbers of final state gluons, we 
will search for a reliable approximation for these exact multiparticle cross sections. 



The Parke- Taylor Formula 

As a first step to such an approximation, we will use as the squared amplitude for 
the process gg — > (n — 2) g the so-called Parke- Taylor (PT) formula ||13|| , which was 



conjectured in Ref. and later proved in Ref. |TT|. Quite surprisingly, the PT result 



can also be derived from soft-gluon factorization techniques ||14|| . This formula is the 
result in leading order in the number of colors (iV c = 3) for the square of the amplitude 
describing the process gg — > (n — 2)g, summed over colors and over a particular set of 
gluon helicity configurations : 

K T r = 2^- 4 ivr 2 (iv c 2 - 1) E4 E — - — a) 

i>j P'(l,2,...,n) S 12S23S 3 4."S„l 

where g s is the QCD coupling constant, = (pi + pj) 2 = 2pi ■ pj and the primed 
sum is over the (n — l)!/2 non-cyclic permutations of (1, 2, n) . The relevant set of 
configurations is the one where all gluons but two have the same helicities, the so-called 
maximally helicity- violating (MHV) amplitudes. There are n{n — 1) MHV amplitudes 
for n > 4. 

This leading N c result contains all the leading collinear and soft divergencies in it 
and furthermore it preserves all the coherence effects and includes initial and final state 
radiation interference. 

In comparing the PT formula with exact results one should devise methods to take 
into account the contributions from amplitudes other than MHV ones. Here we'll adopt 
the simplest approach, due to Kunszt and Stirling [jn| which is to assume that all the 
amplitudes are of comparable size. We'll call this procedure the KS approximation. 
Hence, since there are 2 n — 2(n + 1) non-zero helicity amplitudes in total we simply 
multiply the PT formula by a combinatorial factor which counts the number of non- 
vanishing amplitudes: 

I aKS\2 _ 2" — 2{n + 1) apt\2 / 9 \ 
n(n — 1) 

Several groups have performed comparisons between the PT formula and the results 
from exact calculations. In order to find finite cross sections, a certain set of cuts must 
be imposed on the final state parton momenta. Usually one is interested in large 
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transverse momentum jets which may turn out to be backgrounds for interesting new 
physics. Typical cuts are : pxi > 50 GeV ; \r]i\ < 3 ; 9ij > 30° for the LHC. Kunszt 
and Stirling |T5[ have found good agreement between the exact gg — > gggg cross section 
and the PT formula in the KS approximation and they have also used the PT result 
to compute up to gg — > gggggg and found a broad agreement with the multiplicity of 
large px jets up to six jets seen at the UA1 detector. The exact gg — > ggggg result 



has also been tested against various approximations schemes Hl2| , jig , [17| for different 
sets of cuts. The overall conclusion is that the PT formula in the KS approximation 
overestimates the exact results by as much as 50% for tight cuts but becomes a better 
approximation for looser cuts. This can be understood since for loose cuts the event rate 
is dominated by the infrared and collinear singularities which are taken into account 
by the PT formula. Since all the gluons produced are supposed to be well separated 
in these calculations, the merging of two or more gluons into jets was not taken into 
account. 



Phase Space Parametrization 

The form in which the PT formula is written suggests a parametrization of the (n — 2) 
body massless phase space in terms of the variables s it i+ i = (pi + £>i+i) 2 |[8| : 



(47r)- 2 ™+ 7 1 p dM*_ 3 f4-2, n 



d (PS) = / *=f- / " n-l x 

2 s Jo s - M„_ 3 Jo 



-i 



M u- 3 dM%_ 4 



ML.. _ M 2 J 



'n — 3, n — 2 

"S n _3 n _2 X • • • X 



M l dMl (4a , / 



o Ml -Ml Jo Jo 



ds 34 x / ds 23 (3) 



where Mf = (pi + p 2 H h Pi) 2 , M 2 _ 2 = s , M^_ 1 = M\ = and the limits of 

integration on the variables s» j+i are given by : 



<m = ^l^m - M ?\ ( M ? ~ M li)- ( 4 ) 

There are (n — 3) independent s it i+ i variables, (n — 4) independent Mf variables and 
the integration over the remaining (n — 3) azimuthal angles was already performed. 
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We can write down schematically the cross section for the process gg — > (n — 2)g : 

In order to avoid collinear and infrared divergencies we introduce a Lorentz invariant 
cut-off by requiring Sij = (pi + Pj) 2 > s cut for all particles involved. In particular, 
we will not consider merging processes, i.e., contributions to the n-gluon cross section 
from processes with a larger number of gluons where one or more gluons fail to pass 
our cut. 

We are interested in the large n behavior of this cross section. Due to the factorially 
growing number of permutations entering in the calculation of the PT formula and the 
large dimensionality of phase space, we make two approximations that we feel are 
plausible in the large-n limit : 

• Ordering Independence: For large n, the symmetry of the PT formula and of our 
cuts suggests that each term makes a similar contribution to the cross section 
after integrating over phase space |nj. It is the presence of the two fixed initial 
momenta which violates this hypothesis, and it is plausible that the effect of 
these momenta is less important at large n. In fact, we have verified this approx- 
imation numerically for several arbitrary permutations of the momenta, in the 
cases n = 10, n = 12, using the phase space generators GENBOD and RAMBO. 
This simplifies our calculation tremendously since we now have to compute the 
contribution from the basic string only. 

• Neglect of end point effects: there are (n— 3) independent variables j+i parametriz- 
ing the (n — 2)-body phase space whereas there are n such variables appearing in 

a given string in the PT formula. It is very difficult to write out the 3 dependent 
quantities in term of the independent variables. Since the number of dependent 
variables is n independent, we expect that in the large-n limit their contribution 
is not relevant. Therefore, here we assume that we can conservatively substitute 
s 3 for the dependent variables. 

In addition to these, we also approximate the numerator of the PT formula by (J2i>j s %) ~ 
s 4 . This we call s 4 dominance . 

In the context of these approximations and with the above cuts we find : 

1 1 1 2 n -2(n + l) (4tt)- 2 ™ +7 

(Tgg-+( n -2)g ~ ^ ' 4(^2 - 1)2 ' ( n _ 2 )! ' n(n - 1) 2s 

{2g 2n-4 N n-2 {N 2 _ 1)} . g 4 . {n _ ^ . J_ . J^^/ s) (6) 
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where 



(Ml_ 3 ) mm s - M„ 2 _ 3 i(M2_ 4 ) min M 2 _3 - M„ 2 _ 4 i(Mf) roin Ml - Ml 

s t,-2, n-l ds n -2, n-l f S t-3, n-2 G?S n _ 3i n _ 2 /" S 3 4 rfs 3 4 /" S 2~3 G?S 2 3 



! s n-2, n-l "'Scut s n-3, n-2 -/Scut S3 4 Js C ut S2 3 

with {Mf) min = 2( .^ ! 2) i s cut and {Mf) max = Mf +l - i s cut . We wrote out explicitely 
the different factors entering in Eq. in order to point out their origin. Respectively 
these factors account for : the usual flux factor, the average over helicities and color of 
the initial gluons, the Bose symmetry factor for the (n — 2) identical particles in the 
final state, the KS approximation, the phase-space overall factor, the overall factor in 
the PT formula, the s 4 dominance, the ordering independence approximation and the 
neglecting of end point effects. 

Let us now focus on the estimation of the multiple integral J n _ 2 (A = s cut /s) given 
by Eq. (|7|), which after a trivial integration over the j+i variables can be written as : 

r(xi) ma x d Xl r(x2)max dX2 rfan-4,)max dx n _A , , 1 — X\ . 

I„_ 2 (A) = / / / x ln(- 



(^l)min 1 X\ J (%2)min 1 ^2 J {Xn-4,)min 1 ^n — 4 A 

• H {1 ~ Xl) £~ Xi) ) ■ " ^)(l - *s)) ■ ln(|x 2 (l - x 3 )(l - *<)) ■ ■ ■ 

ln(— x 2 x 3 ■ ■ -x n _ 6 (l - s„_ 5 )(l -x„- 4 )) • ln(— a; 2 x 3 ■ • -a; n _5(l -x„_ 4 )) (8) 
where the new variables Xi are defined by : 

^n-3 M n-4 ^-5 ^2 

Xl = , x 2 = — ^ , x 3 = —2-^ , . . . , x n _ 4 = (9) 

S SX\ SX\X2 SX\X2 ' ' ' Xn—fiXn—5 

and the limits of integration are : 

( \ _ CjKj )max/min „>. 

yXijmax/min l-*- 1 -'/ 
SX\X2 ' ' ' X{—\ 

At this point we could proceed by simply performing a numerical integration of 
Eq. (El) for different values of A = s cut /s; this we did for a number of gluons up to 
n = 14, and the results will be discussed shortly. However, in this section we would 
like to gain some analytic understanding of the cross section for small enough A and 
large n. For these reasons, we explore ways of finding the leading terms of J n _ 2 (A) in 
this limit. 
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Double Leading Log Behavior 

In order to have an approximate analytical form for J n _2(A), we find that the following 
simplifications in the integration limits are useful : 

Xix 2 ■ ■ ■ a*-! 



* 1 * a "-" i - 1 J?L ln m (l - x ) ln n (x) « ln m (l - x) ln"(x) = c(m, n) (12) 

o 1 — x Jo 1 — X 



1 * im2 ~ Xi - 1 — \n m (l-x)n f 1 * ^-\n m (l-x) = —!—\n m+ \l/A) (13) 
o 1 - x v 'Jo 1-x v y m+l v 1 ' K ' 

where the numbers c(m, n) can be easily computed numerically. Within these approxi- 
mations and using MATHEMATICA we are able to find an analytical form for J n _2(A) 
for n up to 11. For n > 11 the computer time required becomes too large. Neverthe- 
less, the analytical approximation provides us with the double leading-logarithm term 
(DLL) in 1/A , i^£(A) for arbitrary n : 



(JL)^ln" n - 7 (l/A) neven 



In-2 (A) = { (14) 
|(i)^ln 2n - 7 (l/A) nodd 

which becomes dominant at small A. The power of In 2 in ([TJ|) reflects the kinematic 
constraint that not every final state gluon momentum can be within s cu t of another 
one. This consideration becomes unimportant at large n. 

We have compared the results from the numerical integration with the approximate 
analytical and DLL results for different values of A and n . We find that the analytical 
and numerical integration results agree over a wide range of n and A. For fixed n, our 
DLL approximation is valid only at small enough A, the range of validity being related 
to the number of partons n. For our results, a rough measure of the range of validity 
of the DLL is given by : 

A < exp[-n] . (15) 

Therefore, from the above considerations we finally arrive at the asymptotic be- 
havior of the small A, large n-gluon exclusive cross section by combining Eqs. (0) and 
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(|HD in the large n limit 



&«\{i$i3"M*r <i6) 

In Eq. ( |T6D the superscript refers to the number of unresolved gluons. We believe this 
result is valid in the range n$ < n ~ n > ln(l/A), where n is the minimum number 
of gluons such that our approximations become reliable. We don't have a definite 
estimate for uq but it is not unreasonable to assume no ~ 10. 



Results for Hard Jets 

Even at moderate values of n and A, where the DLL approximation is not valid, our 
numerical results may be of relevance to SSC physics. This can be seen by writing 
Eq. (|) as : 

a n = - (z(n,A)) n - 2 , (17) 
s 



where 



i 



z(n,A) = -^(I„_ 2 (A))^ . (is) 

When z(n, A) « 1 one would be suspicious of the use of perturbation theory. In Table 1 
we show z(n, A) for some interesting values of n and A, where we conservatively used 
a s (Q 2 = s cut ) = 0.12. Our results suggest that at the SSC, even reasonably "hard" jet 
cuts place A at a value (w 10~ 5 ) which may be approaching a threshold of uncertainty 
for perturbation theory, in the sense that o^+i > o~!£\ In the DLL approximation, 
z(n, A) becomes the effective n-independent expansion parameter of Eq. [Tj| 

.(n,A)^,(A) = -%=ln 2 (l/A) . (19) 

Z7TV iZ 



Discussion and Concluding Remarks 

(1) It is not difficult to understand how our basic result Eq. ([16]) comes about: The 
factor ln 2n (l/A) can be accounted for by simply counting the integrals in Eq. (]?]); it 
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is a consequence of the infrared and collinear divergences contained in the PT formula 
that are regulated by s cut . However, the crucial feature of Eq. fllED is the absence of 
a 1/n! factor. This can be traced to the approximate ordering independence of the 
phase space integrated PT formula; the different n\ permutations compensate for the 
Bose symmetry factor. It should be noted that this cancellation of the n\ does not 
occur for QED: in that case, the n\ terms in the amplitude add to a single term which 



incorporates the permutation symmetry[20] 



(2) The Parke- Taylor formula drops terms in the cross section which are non-leading 



in 1/N C . Non- leading terms in 1/N C do not contribute in DLL, [21|, so that we do not 
consider them further. 

(3) As a result of the absence of the 1/n! or other similar factor, the cross sections 
for the production of n — 2 resolved gluons in the DLL form a geometric series in 

z, and the exclusive total cross section 

"tot — 2.^ °n 
n 

is not summable (except for a phase space cutoff) for z(A) > 1. In particular, there is 
no sign of the exponential summation a ~ exp(— ^ Ins) proposed by Lipatov p2| . We 
will comment shortly on possible implications for inclusive n-jet cross sections. 

(4) The factor y/12 in our small-A result (Eq. |16|) contains the the KS approximation 
(Eq. |]). If this approximation overestimates the cross section by a factor (1 + S) n , then 
the substitution V^12 — > ^12(1 + 5) should be made in Eq. |16[ 

(5) Our results suggest that caution should be exercised before using tree-level exact 
multi-gluon amplitudes to compute mini-jet cross sections. When convoluting the hard 
cross section with the gluon distribution functions, it is possible that (for n not too 
large) the important contributions to the cross section come from low values of the 
initial gluons' invariant mass ^Js {i.e., larger values of A), in which case the use of 
perturbation theory will be valid. However, there will always be potentially large 
contributions from non-perturbative regions of s/s cut , and one must be certain that 
they are under control. 

(6) Can higher order corrections restore the validity of perturbation theory? Schemat- 
ically, we would have to compute contributions to the inclusive cross section o~ n from 
processes with a larger number of gluons where one or more gluons would fail to pass 
our cuts and would be merged to another gluon: 

On = ^ + + 0% + ... (20) 

where the superscripts indicate the number of "unresolved" merged gluons and 
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trf = A <ln 2 "(l/A) 

*#i = A<ln 2 "(l/A)a s (ln 2 (l/e)-ln 2 (l/A)) 

4% = A 2 <ln 2n (l/A)« 2 (aln 4 (l/ e ) + 61n 2 (l/ e )ln 2 (l/A) + cln 4 (l/A)) 

; ; (21) 

e is an infrared cut-off whose contributions will be canceled by loop corrections, and 
the various constants depend on the details of the merging procedure. An explicit 
discussion in the case of e + e~ — > qq + gluons has been given by Brown and Stirling 
and will be commented on below. In our case, 



a n = A a n s ln 2 "(l/A) (l - ^a.ln 2 (l/A) + ^a 2 ln 4 (l/A) + ...)=/„ 

(22) 

In order that perturbation theory be restored with respect to a n , it would be necessary 
that f n z n be bounded for large enough x and n. As an example, the usual n-independent 
Sudakov factor, f n (z) = exp[— Kz\ , is not enough to solve our problem. In their study 
of e + e~ — > jets, Brown and Stirling [23] have examined inclusive jet production using 
a jet merging alogorithm similar to ours (the JADE algorithm), and have not found 
evidence for the presence of Sudakov summation. At this time, we have nothing to say 
about the function f n (z), but its evaluation will be important in proceeding beyond 
the results of this work. 
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A 


n= 10 


n = 11 


n= 12 


ra= 13 


1(T 4 


0.69 


0.70 


0.70 


0.70 


lir 5 


1.1 


1.1 


1.1 


1.1 


10~ 6 


1.5 


1.5 


1.6 


1.7 



Table 1: Numerical values of the effective expansion parameter z(n,A) for different 
values of A and n . 
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